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ABSTRACT 

Objectives: Suicide is a major global health problem 
imposing a considerable burden on populations in 
terms of disability-adjusted life years. There has been 
an increasing trend in fatal and attempted suicide in 
Iran over the past few decades. The aim of the current 
study was to assess overall, gender and social 
inequalities across Iran's provinces during 2006-2010. 
Design: Ecological study. 

Setting: The data on distribution of population at the 
provinces were obtained from the Statistical Centre of 
Iran. The data on the annual number of deaths caused 
by suicide in each province were gathered from the 
Iranian Forensic Medicine Organization. 
IVIetliods: Suicide mortality rate per 100 000 
population was calculated. Human Development Index 
was used as the provinces' social rank. Gini coefficient, 
rate ratio and Kunst and Mackenbach relative index of 
inequality were used to assess overall, gender and 
social inequalities, respectively. Annual percentage 
change was calculated using Joinpoint regression. 
Results: Suicide mortality has slightly increased in 
Iran during 2006-2010. There was a substantial and 
constant overall inequality across the country over the 
study period. Male-to-female rate ratio was 2.34 (95% 
C1 1.45 to 3.79) over the same period. There were 
social inequalities in suicide mortality in favour of 
people in better-off provinces. In addition, there was an 
increasing trend in these social disparities over time, 
although it was not statistically significant. 
Conclusions: We found substantial overall, gender 
and social disparities in the distribution of suicide 
mortality across the provinces in Iran. The findings 
showed that men in the provinces with low 
socioeconomic status are at higher risk of suicide 
mortality. Further analyses are needed to explain these 
disparities. 



INTRODUCTION 

Suicide is considered as one of the three 
leading causes of death among the 15- 
44 years age group and the second cause of 
death in the 15-19 years age group/ It 
imposes a considerable burden on 



Strength and limitations of this study 



■ This is the first national study to evaluate 
regional social inequalities in suicide mortality 
over a 5-year period. 

■ Social inequality in suicide mortality was evalu- 
ated using Cuzick's test for trend and two 
common inequality measures: rate ratio and 
Kunst and Mackenbach relative index of 
inequality. 

■ Age and gender differences between the pro- 
vinces were naively adjusted, which might have 
not fully captured differences in age distribution 
across provinces. 



populations in terms of disability-adjusted 
life years and it has been projected that 
suicide will compose about 2.4% of the 
global burden of diseases by 2020.^ 

Similar to other developing countries, Iran 
has been experiencing a rapid increase in 
suicide rates during recent years. A recent 
study showed that suicide and attempted 
suicide in Iran have increased from 8.3/ 
100 000 population in 2001 to 16.3 in 2007.^ 
Moreover, suicide accounts for 4% of injury 
cases admitted to general hospitals in the 
country.^ 

The risk factors of suicide are some demo- 
graphical characteristics,^ socioeconomic 
situations^ and medical conditions.^ 
However, there are also factors related to the 
area of residence that influence the suicide 
rate.^ Evidence demonstrates persistent geo- 
graphical disparities in distribution of suicide 
between and within countries, which sup- 
ports area-level correlates of suicide.^ 

Iran is a Middle Eastern country with 1 
628 550 km^ land area^^ and consists of 31 
provinces^ ^ that are in different levels of 
development. Ethnic groups tend to reside 
in neighbouring provinces. Therefore, there 
are variations in socioeconomic level and 
ethnicity across provinces as well as 
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geographical and ecological differences that could cause 
a disparity among suicide mortality rates. 

While information on the individual risk factors and 
outcomes of suicide attempts have been well documen- 
ted in previous studies in Iran,^ ^ 12-15 ^Yicyc is little 
information on the role of socioeconomic factors and 
suicide incidence. In particular, there is no known study 
evaluating regional socioeconomic disparities and their 
impact on suicide rates in Iran. 

To fill the gap, this study aimed to describe overall 
and social inequality in suicide mortality rates across all 
provinces in Iran from 2006 to 2010. Although this is an 
ecological study, it will provide a useful starting point for 
examining the social disparity of suicide and provide 
valuable information for policymakers in order to priori- 
tise prevention strategies. 

METHODS 

The data on distribution of population at the provinces 
were obtained from the Statistical Centre of Iran. It 
should be noted that at the time of conducting the 
current study, Iran constituted of 30 provinces; it was 
later that Tehran province was split into two provinces. 
The data on the annual number of deaths caused by 
suicide in each province were gathered from the pub- 
lished reports of Iranian Forensic Medicine 
Organization, which is affiliated to the Judicial Authority 
in Iran. According to Iranian law, all deaths due to exter- 
nal causes should be reported to this organisation for 
examination and recording and to issue the death cer- 
tificate. It is considered as the most reliable source of 
mortality data in Iran.^'^ Next, annual suicide mortality 
rates per 100 000 population were calculated for each 
province. Human Development Index (HDI) was used 
as the provinces' social rank and related data were 
obtained from the President Deputy of Strategic 
Planning and Control. The HDI is a composite index of 
three basic dimensions of human development, includ- 
ing life expectancy at birth, educational attainment 
(based on a combination of adult literacy rate and 
primary education to tertiary education enrolment 
rates) and income (based on gross domestic product per 
head adjusted for purchasing-power parity (US$))^^. As 
a composite index, it is expected that HDI might 
capture socioeconomic development more comprehen- 
sively than single indicators such as average income or 
expenditures. 

Overall inequality measures inequalities in health, irre- 
spective of the other characteristics of the individual. 
To measure overall inequality we followed the same 
approach as measuring income inequality and used the 
Lorenz curve and Gini coefficient. These two measures 
are commonly used in assessing overall inequality in dis- 
tribution of healthcare resources and outcomes. 
The Lorenz curve is used to compare the distribution of 
health measure with perfect equality (diagonal line). In 
the current study, the Lorenz curve was plotted as the 



cumulative share of population ranked by suicide mor- 
tality rate, in an increasing order, against the cumulative 
share of suicide mortality. The further the distance from 
the diagonal line the greater the degree of inequality. 
The Gini coefficient is equal to twice the area between 
the Lorenz curve and the diagonal line and its value 
ranges from 0 (perfect equality) to 1 (maximum pos- 
sible inequality). This measure takes into account the 
distribution of health variables across the entire popula- 
tion. In the current study, we used fastgini command in 
STATA to calculate the Gini coefficient and its jackknife 
95% CI. In order to examine the gendered nature of 
suicide mortalities (ie, gender inequality) in Iran, we cal- 
culated the male-to-female rate ratio (RR) and its 95% 
CI using negative binomial regression with a robust vari- 
ance. As we only had data stratified by sex groups for the 
whole study period and not for every specific year, an 
overall RR was calculated. 

Social inequality was evaluated using Cuzick's test for 
trend and two common inequality measures: RR and 
Kunst and Mackenbach relative index of inequality 
(RIIkm)-^^ Cuzick's test for trend is an extension of the 
Wilcoxon rank-sum test and is used as a non-parametric 
test for trend across ordered groups. To calculate RR, 
the provinces were ranked and divided in five quintiles 
by HDI (weighted by their population). Then, negative 
binomial regression with a robust variance was used to 
calculate RR and its 95% CI to compare the highest 
versus the lowest quintile. One problem with RR is that 
it only considers the population in two extreme socio- 
economic groups. To take into account the whole popu- 
lation distribution across socioeconomic groups and also 
to remove differences in the size of socioeconomic 
groups, as a source of variation in the magnitude of 
health inequalities, RIIrm was calculated.^^ RII is widely 
used to measure social inequality and is recommended 
when making comparisons over time or across popula- 
tions.^^ To calculate RII, after determining the relative 
position of the population in the provinces ranked by 
HDI, the number of deaths in the provinces was 
regressed on these relative ranks using negative binomial 
regression with a robust variance and population as 
offset variable. With the lowest social rank as reference, 
an RIIkm value greater (lesser) than 1 shows that the 
suicide mortality rate was higher among the provinces 
with higher (lower) social rank (more distance from 1 
implies more inequality) .^^ To account for sex and age 
differences between the provinces, we also estimated 
adjusted RIIrm by including the proportion of males in 
the population, and the mean age of females and mean 
age of males in our negative binomial regression. 

To examine temporal changes of suicide mortality 
rates across the provinces and also across the quintiles of 
HDI, we calculated annual percentage change (APC) 
and its 95% CI for each province and quintile using the 
Joinpoint Regression Program V.3.5.4. Moreover, this 
program was used to calculate APC of overall and socio- 
economic inequality measures over the study period. 
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APC is estimated using the following regression: 

Ln(It) =bo+bi(t) 

APC = (e^i - 1) X 100 
where It shows the suicide mortality rate and estimated 
inequality indices for year t. 

In a sensitivity analysis, Tehran was excluded from the 
analysis to examine the overall, gender and social 
inequalities across the remaining provinces. The reason 
for this was that Tehran has special status as the capital 
of the country and is the centre of economic, social and 
political activities. Excel office and STATA V.ll 
(StataCorp LP, College Station, Texas, USA) were used 
for statistical analysis. 

RESULTS 

Table 1 shows mean population, HDI, suicide mortality 
rate and APC (%) across the provinces for years 
2006-2010. Ham and Hormozgan provinces had the 



highest and the lowest suicide mortality rate during the 
study period, respectively (8.8-fold difference). Most pro- 
vinces (56.6%) had a suicide mortality rate of 3-6/100 000 
population. In addition, suicide mortality rate was more 
prevalent among the western provinces of Iran (figure 1). 

The estimated APC values show that only five pro- 
vinces experienced significant changes in suicide mortal- 
ity rate over the study period (significant increases in 
Ham, Isfahan, North Khorasan and Tehran provinces 
and significant decrease in Markazi province). 

Figure 2A shows temporal changes of suicide mortality 
rates for the country. While the graph shows a slight 
increase in suicide mortality rates over the study period 
(from 4.25 in 2006 to 4.88 in 2010), this was not statistic- 
ally significant (APC=3.05%, p=0.23). Figure 2B presents 
suicide mortality rates for five quintiles of HDI and 
examining the trends showed that APC was significant 
only in the highest quintile (APC=17.98, p=0.01). 
Figure 3 shows scatter graphs of HDI and suicide mortal- 
ity rates. It is evident that the higher HDI was associated 



Table 1 Mean population, Human Development Index (HDI), suicide mortality rate and annual percentage change (APC) 
across the provinces in Iran, 2006-2010 (ranked by suicide mortality rate) 

Suicide mortality 

I rate per 100 000 





Population 


HDI 


population 


APC (%)* 


Nam 


555 929 


0.729 


19.53 


10.11 


Kermanshah 


1 892 100 


0.748 


13.74 


0.46 


Lorestan 


1 736 946 


0.761 


10.64 


1.18 


Hamedan 


1 700 960 


0.740 


9.59 


-1.52 


Gilan 


2 428 553 


0.769 


6.52 


9.98 


Ardebil 


1 235 234 


0.735 


6.22 


-2.32 


East Azerbaijan 


3 646 459 


0.763 


5.58 


1.31 


Zanjan 


973 739 


0.752 


5.53 


12.14 


Kohgyluyeh and Boyerahmad 


651 577 


0.718 


5.25 


-4.20 


Khuzestan 


4 372 242 


0.762 


5.11 


-7.13 


West Azerbaijan 


2 944 224 


0.713 


4.72 


-2.74 


Kordestan 


1 453 503 


0.713 


4.69 


5.25 


Golestan 


1 651 708 


0.737 


4.49 


5.67 


Overall (Iran) 


72 599 045 


0.758 


4.46 


3.05 


Qazvin 


1 177 582 


0.783 


4.25 


-1.13 


North Khorasan 


824 979 


0.759 


4.24 


10.64 


Pars 


4 431 684 


0.783 


4.16 


1.56 


Chaharmahal Bakhtiari 


875 207 


0.749 


4.00 


-7.72 


Qom 


1 087 01 1 


0.773 


3.75 


7.89 


Mazandaran 


2 979 189 


0.745 


3.66 


-5.71 


Isfahan 


4 680 831 


0.810 


3.64 


4.62 


Kerman 


2 799 417 


0.750 


3.10 


-5.03 


Bushehr 


914710 


0.786 


3.02 


15.55 


South Khorasan 


656 469 


0.723 


3.02 


-2.28 


Tehran 


14 106 297 


0.843 


2.92 


17.98 


Yazd 


1 028 152 


0.809 


2.68 


-4.48 


Razavi Khorasan 


5 765 706 


0.777 


2.66 


-2.64 


Semnan 


606 982 


0.814 


2.60 


4.27 


Markazi 


1 371 514 


0.785 


2.38 


-10.73 


Sistan and Baluchestan 


2 569107 


0.643 


2.23 


6.11 


Hormozgan 


1 481 031 


0.766 


2.21 


-5.59 



*ltalic figures show statistically significant results (p<0.05). 
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Figure 1 Distribution of suicide mortality rates across the provinces in Iran. 



Male-to-female RR was 2.34 (95% CI 1.45 to 3.79), 
implying a significant higher suicide mortality rate 
among males than females over the study period. 
Examining this ratio in the provinces showed that in all 
provinces but Kordestan men had a higher suicide mor- 
tality rate than women. Excluding Tehran from the 
sample did not change this finding (2.33; 95% CI 1.43 
to 3.79). 

Table 2 presents the overall and social inequality mea- 
sures in distribution of suicide mortality rates across the 
country through 2006-2010. The Gini coefficient 
ranged from 0.248 to 0.302, implying substantial overall 
inequality across the provinces. The Lorenz curves corre- 
sponding to this Gini coefficient is shown in figure 4. 
There was a 14.5% decrease in the Gini coefficient 
between the first year and the last year of study, implying 
decreasing overall inequality between these two points of 
time. Over the study period, the APC of the Gini coeffi- 
cient was —4.28 and statistically non-significant. 
Excluding Tehran resulted in a 6.9% increase in the 
Gini coefficient between the first and the last year of the 
study period and a statistically non-significant APC of 
1.42 was estimated. 

RR was significantly lower than 1 in overall and for all 
years of the study period, implying a higher suicide mor- 
tality rate among people in the lowest quintile of HDI 
compared with the highest. Over the study period, RR 
was approaching 1 and the APC of RR was 12.2 and 



with lower suicide mortality rates and Cuzick's test con- 
firmed this (z=-4.61, p<0.011). 
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Figure 2 Suicide mortality rates per 100 000 populations in 
(A) total sample; (B) quintiles (Q) of Human Development 
Index, 2006-2010. 
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Figure 3 Scatter plots of Human Development Index and suicide mortality rates, stratified by year of study. 



statistically significant, implying a decrease in the gap 
between the highest and the lowest social ranks. 
Although excluding Tehran from the sample did not 
change the overall picture of social inequality, an inverse 
trend (increasing social inequality) was observed. 

RIIkm was also lower than 1 in overall and for all years 
of the study, showing a persistent inequality in favour of 
people living in the provinces with higher social rank. 
Temporal analysis showed that RIIrm did not signifi- 
cantly change over the study period. Adjusting for age 
and sex did not change this observation. When we 
excluded Tehran province, the APC value for adjusted 
RIIkm was statistically significant, showing an increase in 
social disparity across the remaining provinces in Iran. 

DISCUSSION 

In this study, for the first time, we assessed overall and 
social disparities in the distribution of suicide mortality 
across the provinces of Iran over a period of 5 years 
(2006-2010). The findings showed that suicide mortality 
has slightly increased over the study period. The findings 
also indicated that there were substantial overall, gender 
and social disparities in the distribution of suicide mortal- 
ity across the country and it was higher in the provinces 
with lower social rank. These disparities were generally 
stable and persistent over the study period. 

The findings from the current study showed an 
inverse association between the provinces' social rank 
and suicide mortality in Iran. Although the studies on 
association between area social rank and suicide mortal- 
ity reported mixed results,^ the results of this study 



are in line with previous ecological studies investigating 
the relationship between socioeconomic characteristics 
and suicide rates, in particular those studies focused on 

Q ofy go Q 

high-income settings. Rehkopf and Buka, in their 

systematic review of suicide and socioeconomic 
characteristics of geographical areas found that, among 
studies with statistically significant results, 50% and 73% 
of studies reported an inverse relation between area 
income and education characteristics and completed 
suicide, respectively. They also found that the probability 
of reporting an inverse relationship between area social 
rank and suicide mortality was higher among the studies 
conducted in Asia (94% of studies with statistically sig- 
nificant results).^ Similar findings were reported by 
another study^"^ focusing on countries in the Eastern 
Mediterranean region (where Iran is located) , indicating 
that high-income countries had lower suicide mortality 
rates than their low-income and middle-income counter- 
parts. It is argued that people in the provinces with 
lower social rank generally have more adverse experi- 
ences, poorer mental health, lower access to psychiatric 
services and lower access to health facilities. These 
factors might partly explain higher suicide mortality rate 
in the provinces with lower social rank in Iran. 

The four Western provinces of Iran (ie. Ham, 
Kermanshah, Lorestan and Hamedan) had the highest 
suicide mortality rate in the country. One potential 
explanation for this observation can be the low socio- 
economic status of these provinces. These provinces are 
among the provinces with the highest unemployment 
rate and lowest HDI in the country. High divorce rates 
in these provinces (except for Ham) can be another 
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potential explanation, as divorce is considered as a risk 
factor for suicide mortality. In addition, cultural 
issues such as the tribal structure of communities and 
the extreme fanaticism prevailing in these provinces 
have been considered as another potential explanation 
for this finding.^^ 

The high gender gap in suicide mortality rates 
observed in the current study is in line with the previous 
epidemiological studies in Iran^ ^ and is comparable to 
the studies conducted in other settings, in particular 
high-income countries. Although many studies, 

including the previous studies in Iran,^ ^ have shown 
higher suicide attempt rates among females, completed 
(fatal) suicides are higher among males. One potential 
explanation can be the difference in methods of attempt- 
ing suicide among males and females. For example, the 
most common methods of attempting suicide used by 
males in Iran are by hanging and using firearms, which 
have higher fatality rates compared with the self-burning 
method commonly used by females.^ Greater psy- 

chosocial impact of problems, such as unemployment or 
retirement, on males compared with females^^ and 
adopting coping strategies such as emotional inexpres- 
siveness, lack of help-seeking, risk-taking behaviour, vio- 
lence, alcohol abuse and drug misuse by males (which 
are triggered by norms of traditional masculinity) are 
other potential factors that have been discussed in the lit- 
erature. The findings of the gender analysis are import- 
ant for designing and implementing suicide prevention 
strategies as the factors, patterns and behaviours asso- 
ciated with suicide are affected by gender. 

Temporal analysis of suicide mortality showed interest- 
ing results, indicating that among five quintiles of HDI it 
was only the highest quintile that experienced a signifi- 
cant change in suicide mortality over the study period. 
This finding can potentially be explained by the increas- 
ing prevalence of mental disorders, rising unemploy- 
ment rates and an increasing trend in divorce 
rates^^ in the provinces with higher HDI, such as 
Tehran, the capital city of Iran. 

To our knowledge, this is the first national study evalu- 
ating social inequalities across different regions over 
time. Although this study was conducted in the context of 
Iran, the findings may also be applicable to other 
middle-income countries, in particular countries in the 
Middle East region, which share a similar culture. 
Moreover, we believe that, in terms of methodology, our 
analysis presents a good example for employing a triangu- 
lation of different methods for evaluating inequalities in 
suicide mortalities. However, the current study has several 
limitations that should be considered when interpreting 
its findings. First, age and gender differences between 
the provinces were naively adjusted, which might not 
have fully captured differences in age distribution across 
provinces. Second, there is the issue of availability and 
quality of data on suicides, which is common in develop- 
ing country settings^^ and leads to underestimation 
and misclassification of suicide. This might be an issue in 
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Figure 4 Lorenz curves of the distribution of suicide 
mortality in Iran, 2006-2010. 

this study because of the social stigma of suicide and reli- 
gious sanctions and some legal issues in the Iranian 
context.^ Moreover, the underestimation and misclassi- 
fication of suicide mortality might be more common in 
the provinces with lower social rank; therefore, we expect 
the social disparity to be more profound than what has 
been reported here. Third, the current study is an eco- 
logical study using province as unit of analysis, which 
embraces substantial heterogeneity within provinces. 
This implies that the observed disparities in suicide mor- 
tality are between provinces and they are not necessarily 
applicable to smaller geographic units or individuals. No 
causal inference can be drawn due to the ecological 
nature of the study and, furthermore, there was no 
control for confounders in this study. 

CONCLUSION 

The present study indicated that there were substantial 
overall, gender and social disparities in the distribution 
of suicide mortalities across the provinces in Iran. 
Moreover, the study showed an inverse association 
between the provinces' social rank and suicide mortality. 
The findings imply that prevention resources should be 
targeted in high-risk groups, in particular men in the pro- 
vinces with low socioeconomic status. Further investiga- 
tions are needed to explain these disparities in suicide 
mortality across the provinces. Moreover, more studies 
are needed to explore the association of socioeconomic 
factors and suicide (attempted and fatal), focusing on 
smaller geographical units and at the individual level, in 
order to design better prevention strategies. 
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